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INFINITE PRODUCTS OF ANALYTIC MATRICES 



GEORGE D. BIRKHOFF 

In a large part of the theory of functions of a single complex variable the 
matrix of analytic functions rather than the single analytic function must be 
taken as the fundamental element. This is certainly the case for the functions 
defined by linear difference and differential equations. 

The goal of the present paper is to show that the classical results of Weier- 
strass and Mittag-LefHer, treating of the formation of infinite products of 
functions with assigned singularities, admits of a natural extension to infinite 
products of matrices. The matrices considered will be square matrices of 
rfi elements and of determinant not identically zero. 

The concept of equivalence, which I have developed elsewhere, f lies at 
the basis of this extension: Let 

A {x) = \aij {x)\ and B{x) = {hij{x)\ (i,j = i,---,n), 

be two matrices of analytic functions, each analytic in the neighborhood of 
a point X = Xo but not necessarily analytic at the point Xq. If the matrix 
M (x) defined by the matrix equation 

Aix) = M{x)B{x) 

is composed of elements m^- ( a; ) each analytic aXx = xq and if the determinant 
I If ( X ) I of the matrix M (x) is not zero at a-o , then A{x) is equivalent to 
B {x) at a; = xo (from the left).t 

From this definition it follows immediately that if ^ ( x ) is equivalent to 
5 ( a; ) at a; = Xo , then B (x) is also equivalent to A{x); if, further, B (x) 
is equivalent to C (x) , then A{x) is equivalent to C ( x ) . 

Two matrices equivalent to each other at a point have essentially the same 
type of singularity at the point. It is this fact which gives importance to the 
notion. 

Evidently any two matrices of functions analytic at xo are equivalent if the 
determinants of the two matrices do not vanish at the point. For this reason 
a point for which the elements of a matrix are analytic, and the determinant 

* Presented to the Society, December 28, 1915. 

t See, for example. Proceedings of the American Academy of Arts 
and Sciences, vol. 49 (1913), pp. 521-568; in particular p. 540. 

t For the elements of the theory of matrices assumed in the present paper see Schlesinger, 
Vorlesungen liber linearen Differenlialgleichungen, pp. 18-19. 
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of the matrix is not zero, will be called an ordinary point of the matrix. If 
the elements of the matrix are analytic, while the determinant vanishes, the 
matrix will be said to have a singular point of analytic type at x = a;o; if the 
determinant vanishes to the kth order, this singular point will be said to be 
of multiplicity k . 

If the elements have poles of the A:th order at most at a; = a;o , the matrix 
will be said to have a pole of the kth order; at such a point the determinant of 
the matrix will be analytic or have a pole of order not greater than nk . 

Finally, if the elements are analytic in the neighborhood of a given point 
while some elements have an essential singularity at the point, and if the 
determinant of the matrix does not vanish in the neighborhood of the point, 
the matrix will be said to have an isolated essentially singular point there; if, 
however, the determinant does vanish in every neighborhood of the point, 
or if the point under consideration is a cluster point of singular points of the 
matrix not of analytic type, the matrix will be said to have a non-isolated 
essentially singular point there. 

It is to be recalled that in the theory of the factorization of an analytic 
function a cluster point of zeros does not play the role of an isolated essentially 
singular point. 



1. Factorization of polynomial matrices 

We will commence with the consideration of matrices whose elements are 
polynomials. By the order of such a matrix we shall mean the degree of its 
polynomial determinant in x. 

From the standpoint of factorization the simplest type is that of order zero 
(constant determinant). The inverse of such a matrix and the product of 
two such are also of order zero. These matrices are not singular in the finite 
plane.* 

A matrix of polynomials will have a single singular point of analytic type 
in the finite plane if its determinant is of the form c{x — .To)* . This point 
will lie at X = xo and be of multiplicity k . 

An elementary matrix will be a matrix of the form 



(1) 



U{x) = 



rl 





Ci 


0-1 





1 


Ci 











■ ■ ( X — Xo ) • • 





-0 





Cn 


1- 



* Questions of factorization for such matrices have been considered by Hensel and Lands- 
berg, Theorie der algebraischen Funklionen einer Variabeln, etc., p. 163. 
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and differs from the unit matrix only in some ^-th column where the successive 
elements are the following: 



Cl, Ci, 



, Ck-l, X — Xq, Cjt+1, 



Its only finite singular point will be one of analytic type and multiplicity 1 at 
X = Xq. The kth column of the matrix U (x) will be termed its character- 
istic column. 

Lemma. The most general matrix F (x) with a singular point of analytic 
type at x = Xo of multiplicity k may he decomposed into the product of a matrix 
Q{x) , for which Xo is an ordinary point if k = I or a singular point of analytic 
type and of multiplicity k ~ 1 if k > 1 , and a second elementary matrix U (x) 
with singular point at Xq , i. e., 



(2) 



Fix) = Q{x)U{x). 



Proof.- — The determinant of F (xq) is zero by hypothesis and thus we may 
determine a matrix of constants D of determinant not zero so that the product 
matrix F {xo)D has the elements of some Ith column all zero. To effect this 
it is necessary merely to determine the elements du, d2i, • • • , dni (not all zero) 
of that column of D so that the n equations 

fix (Xo) dii + fi2 (Xo) d2l + • • • + fin iX(,)dnl = 



(i = l,2, 



,n) 



hold, and then to choose the remaining elements of D so that its deteminant 
does not vanish. But this set of n linear homogeneous equations in du, 
dii, • • • , dni can always be satisfied since the determinant of the coefficients 
is |i^(a;o)| = 0. 

Furthermore let us assume for the moment that the first of these quantities 
du is not zero and let us take 1 = 1. In this case the elements of D occurring 
in the n — \ columns after the first may be taken to be 1 along the principal 
diagonal and zero elsewhere. Also the elements of the first column may be 
taken to h& 1, d^, dz, • • • , d„ so that 



D = 



1 

d2 1 



[dn 



1 



When D has been taken in this way, the matrix F {x)D has the property that 
the elements of the first column vanish at x = Xq. In consequence, if we 
write 

X — a-o • • • 

1 ••• 



F{x)D = q{x) 











1 
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SO that Qix) is derived from F {x)D by dividing the elements of the first 
column of F {x)D by x — xq, the matrix Q{x) is made up of elements 
analytic at Xq . 

Moreover, li k = 1 , the determinant of Qix) is not zero at x = a-o; this 
follows from the determinant relation |i^(a;)| = \Q{x)\(x — Xq) . Likewise, 
if fc > 1 , the same relation shows that Q(x) has a singular point of analytic 
type of multiplicity k — 1 . 

If we multiply the above matrix equation through by D~^ on the right, we 
obtain 

X — Xo • • • 

Fix) = Qix) -^ 1 ■•• 



- dn • • • , 



The second factor is evidently an elementary matrix Z7(x) while the first 
factor Qix) has the properties demanded of it in the lemma. The first 
column of U ix) is the characteristic column. 

In the case when du is zero some one of the quantities dii, • • • , dni will be 
different from zero. If d^i is not zero, for example, we may take / = 2 and 
factor F ix) into the product of a matrix Qix) and an elementary matrix 
U ix) whose second column is its characteristic column. In general if the 
Xth of the quantities du, • • • , dni is not zero we may take I = X and proceed 
as before. 

This completes our proof of the lemma. 

It is evident, in the simple case when Fix) has a singular point of multi- 
plicity one at a; = a-o , that there may be n distinct ways of separating out an 
elementary factor, namely by allowing / to range from 1 to n. Hence the 
theory of factorization for w. > 1 diverges greatly from that for the case 
n = 1 inasmuch as the mode of factorization is not unique. 

The matrix Qix) of the lemma has the same singular points in the finite 
plane for a; + .To as P ( a; ) . Furthermore such a singularity is of the same 
kind for both matrices, i. e., is of analytic type with the same multiplicity, or 
a pole of order k , or an isolated essential singularity, or a non-isolated essential 
singularity for both. This is an obvious consequence of (2) and the inverse 
relation Qix) = Fix)U~^x) . 

We are now in a position to establish at once the following theorem: 

Theorem I. An arbitrary matrix of 'polynomials Fix) of order m whose 
determinant is factorable into the product 

Cix — Xi) ix — X2) ■ ■ ■ ix — Xm) ( C + ) 

may be expressed as a product 
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where Po(x) is a polynomial matrix of order zero and where Ui{x) , TJ2{x) , 

■ ■ ■ , Urn ( X ) are elementary matrices with singular points at x\, Xi, ■ ■ ■ , Xm 
respectively. 

It is obvious that a successive application of the above lemma enables us 
to write P (x) in the stated form. 

Suppose now that a polynomial matrix P (x) be expressed in any way as a 
product Fi ( a- ) F2 ( a- ) • • • Vm{x) of m polynomial matrices Vi{x) , F2 ( a; ) , 
• • • , Vm (x) , each having its only finite singularity at distinct points Xi, X2, 

■ • • , Xm respectively. A particular product of this sort can be obtained at 
once from the product for P (x) of Theorem I. 

The most general factored form Vi(x)Vi{x) • • • F^ ( x ) of P (x) with the 
same properties is defined by the equations 

V^{x) = V^{x)R,{x), 

V,{x) = Rr'{x)V,{x)R2{x), 

F™(x) = Rzliix) Vm{x) , 

where Ri{x) , • • • , iC-i (x) are arbitrary polynomial matrices of order zero. 

It is clear that the indicated general factored form is a possible form. 

Further, any factorization Fi ( a; ) F2 ( a; ) ■ • • Vm{x) leads at once to the 
matrix identity 

V,{x) ■■■ F^Cx) ?„-;>( a;) ••• n' (x) = Ff^a;) Fi(a;) . 

The matrix Ri{x) , as given by the right-hand member, has its only possible 
finite singularity at xi, while, as given by the left-hand member, it has its 
only possible singularities at X2 , • • • ,Xm- Hence i?i ( x ) is not singular in the 
finite plane and must be composed of polynomial elements. This gives us 
the first equation. 

By use of the relation just obtained we find 

F2(x) ••• Vm(x) = R^{x)V,(x) ••• F„(x), 

and from this we get further 

F3(x) ••• V,n{x)V-'{x) ■■■ V^'ix) = F2-'(a;)i?i(x)F2(x). 

The matrix given by either side of this equation must, for a similar reason, 
be a polynomial matrix Ri{x) of order zero. This conclusion leads us at 
once to the second of the given equations. 

Continuing in this way we derive all of the desired equations in succession. 

It is interesting to consider a particular resolution in elementary factors 
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possible in general when the elements of P (x) are polynomials of the juth 
degree and when | P ( a; ) | is a polynomial of degree rifx = m with distinct 
roots a-i , X2 , • • • , Xm- 

If we remove an elementary factor corresponding to x = Xm for which the 
first column is the characteristic column, we thereby only alter the first column 
of P ( a; ) . In fact we obtain the quotient Q,{x) of the lemma from P (x) 
by multiplying on the right by U~^ (x) , which is equivalent to adding a con- 
stant multiple of the second, third, • • •, last column to the first column, and 
dividing that column by x — Xm- Thus the degree of the first column is 
lowered by unity, while the other columns are unaltered. 

In the same way we may lower the degree of the second column by unity, 
by the removal of an elementary factor corresponding to a; = Xm-i for which 
the second column is the characteristic column. 

By this process, after n steps, we lower the degree of all n columns of 
P ( X ) by unity by removing n elementary factors. 

Repeating the process we can reduce the quotient to a constant matrix in 
tin steps. Every step can be carried out in the general case, i. e., if certain 
particular algebraic relations between the coefficients of P ( x ) do not hold. 

Thus the matrix Po ( x ) of Theorem I of order zero can be taken to be a 
matrix of constants in this case. 

This result shows that one can not expect to obtain a simpler normal form 
for a factor with one singular point than the elementary matrix. For we have 
used the tifj, roots Xi, • • • , Xm, the nix elementary factors each involving n — 1 
arbitrary constants, and a constant matrix with n^ constants. Thus there are 
involved ra^ (ju + 1) explicit constants in the factorization of P (x) , precisely 
as many as there are in the polynomial matrix P ( x ) . No further simpli- 
fication can be expected since no superfluous constants have been employed. 

2. Equivalence problem for polynomial matrices 

In the first paragraph we solved the simple problem of factorization of 
polynomial matrices. The converse problem of constructing a polynomial 
matrix with given factors is solved merely by writing down the product of the 
factors. 

But it is possible to modify the point of view slightly. Every polynomial 
matrix with singularities at the distinct points Xi , • • • , x^ is by the lemma of 
§ 1 equivalent to a polynomial matrix Vi ( x ) at the point Xi whose only finite 
singularity lies at Xi for i = 1 , 2 , • • • , m . 

We are thus led to formulate a converse problem which may be stated as 
follows : 

Equivalence problem. Give7i polynomial matrices Fi ( x ) , • • • , F„, ( .r ) , 
each having its only finite singular point at the distinct points xi , • • • , x^ re- 
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spedively; to construct a polynomial matrix P (x) equivalent to these matrices 
atxi, • • • ,Xm respectively hut with no finite singular points other than X\, ■ • • ,Xm- 

Before proceeding to its solution we may note that the formulation of the 
problem is symmetrical with respect to all the singular points involved. 

An explicit solution may be obtained by the aid of a simple notation. Let 
\T {x)]^ denote a polynomial matrix, with its only finite singularity at ^, 
to which the polynomial matrix T {x) is equivalent at the singular point 
a; = ^ofr(a;). The lemma of factorization of § 1 shows that, since T {x) 
is composed of elements analytic at a; = ^ , a matrix [T {x)]^ given as a product 
of elementary matrices, will always exist. If now we write for brevity 



(3) 



Pi{x) = [Fi(.r)U= Fi(x), 
P.ix) = {V,{x)P^'{x)],,, 



Pn.{x) = [V^{x)Pt'{x) ■•■P-U{x)],^, 
a solution is 
(4) P(a;) =P„(x)P„_i(x) •••Pi(x). 

To prove that P {x) , as defined in this manner, has the desired properties 
we begin by observing that P ( x ) is certainly equivalent to Fi ( x ) at x = Xi 
since the product Pm ( x ) • • • P2 ( x ) is not singular at x = xi . 

Next let us proceed to the point X2 . By definition of P2 ( x ) we may write 

P2(x) = ^(x)F2(x)Pr'(x) 

where A{x) is not singular at x = X2 . Consequently, on substituting in 
this value for P2 ( x ) , we see that 

P(x) =P„(x) •••P3(x)^(x)F2(x). 

But Pm ( X ) , • • • , P3 ( X ) and A{x) are not singular at x = X2 so that P ( x ) 
is evidently equivalent to Vi ( x ) at x = Xi. 

In precisely the same way we proceed with the proof that P ( x ) is equiv- 
alent to F3 ( X ) , • • • , Fm ( X ) at X3 , • • • , Xm respectively, and thus forms a 
solution. 

The most general solution P ( x ) of the equivalence problem is Po ( x ) P ( x ) 
where Po(x) is any polynomial matrix of order zero. 

In the first place it is clear that P ( x ) as so defined is a solution of the equiv- 
alence problem. Moreover, consider any solution P ( x ) . At any singular 
point as Xi we have 

P(x) = ^(x)Fi(x), 

P(x) = P(x)Fi(x), 

where A{x) and B{x) are not singular at .t = xi . Using these values of 
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P and P we conclude that the matrix P ( x ) [ P ( a* ) ]~^ is equal to B{x) A~^ {x) 
and thus is not singular at Xi. Similarly this matrix is not singular at x^, 

• • • , Xm and also of course at any other point of the finite a--plane. Since 
P ( X ) by hypothesis is a polynomial matrix the matrix Po ( x ) under consider- 
ation must be a polynomial matrix of order zero, which we were to prove. 

3. Factorization of matrices of entire functions. First method 

Let E{x) be a matrix of entire functions. Its finite singularities are all of 
analytic type, and fall at the points Xi, .X2 , • • • , at which the determinant of 
E{x) vanishes. We take these points to be arranged in some arbitrary order 
and assume that they are infinitely many in number. It is clear that we have 
lim Xm = » . 

Suppose that Fi (x) is some polynomial matrix with unique finite singularity 
at Xi , to which E (x) is equivalent at xi; for instance Vi{x) may be a product 
of elementary matrices with singularity at Xi , as follows from the lemma of 
§ 1. The matrix E{x) Ff' (x) is a matrix of entire functions without a sin- 
gularity at X = Xi but with its other finite singularities at X2 , X3 , • • • . At X2 
this modified matrix is equivalent to a polynomial matrix F2 ( x ) with unique 
finite singularity at X2 . Forming now the matrix of entire functions 
E(x) F7' ( X ) Fj'' ( X ) , which is not singular at xi and X2 , we are led to F3 ( x ) , 
V4ix) , • • • in similar fashion. 

It is therefore possible to find such polynomial matrices Fi ( x ) , 1^2 ( x ) , 

• • • that E{x) is equivalent to Fi ( x ) at Xi , to F2 ( x ) P''i ( x ) at X2 , and so on. 
The following theorem shows that matrices E{x) of entire functions having 
this property exist for any choice of Fi ( x ) , F2 ( x ) , • • • with a single finite 
singularity of analytic type at Xi , X2 , • • • respectively. 

Theorem II. If Vi (x) , F2 (x) , • • • are polynomial matrices each ivith a 
single finite singularity of analytic type at the distinct points xi , X2 , • • • and if 
lim Xm — 00 , there exists a matrix of entire functions equivalent to the product 
Vm{x) • ■ ■ Fi ( X ) at Xm (m = 1,2, • • • ) and of the form 

(5) E (x) = lim Qm{x)Vm{x) • • • Fi (x) , 

m=oo 

where Qm{x) (m = 1, 2, •■ ■) is a matrix of entire functions without singxi- 
larities in the finite plane. 

Proof. Let us first attempt to choose Qi(x), Qiix), ••• so that the 
sequence of matrices 

Pm{x) = Qmix)Vm{x) • • • Fi ( X ) 

converges to a limit matrix as m approaches infinity. Now we have clearly 

(6) Rm+iix) = [I + Smix)][I + Sm-iix)] ■■■ [ 7 + Si ( X ) ] Qi ( X ) , 
where I + Sm{x) = Q^+i ( x ) Vm+i ( x ) Q''^ (x) . 
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It will therefore suffice to establish that the matrices Sm(.x) can be made 
to have elements of as small absolute value as desired within a circle of radius 
less than Xm with center at the origin. This will be established by means of 
the following lemma: 

Lemma. If F (x) be a matrix without singularity in some closed continuum C 
in the finite plane, it is iwssihle to find a matrix Fi{x) without singularity in the 
finite plane and such that the absolute values of the elements of the difference 
matrix F (x) — Fi{x) are arbitrarily small throughout C . 

A very simple proof of the lemma may be given as follows. The elements 
of the matrix 

dF{x) 



dx 



-F-Hx) = A(x) 



are not singular in C . The elements ay ( x ) of ^ ( a; ) may then be represented 
by series of polynomials which converge uniformly throughout C , by a well- 
known theorem due to Runge.* Hence we may break off the terms of these 
series at such a place that the elements of the resulting matrix Ai{x) differ 
from those of ^ ( cc ) by a uniformly small quantity throughout C . 
Consider now the matrix differential equation 

dFij x) 

— ^— = Aiix)Fiix). 

By comparison with the equation for F , 

dF{x) 



dx 



= A(x)F{x), 



we infer that a solution Fi{x) of the preceding equation may be found whose 
elements are analytic and differ from those of F (x) by as little as we please 
throughout the given continuum C . 

But the differential equation for Fi is without finite singular points so that 
|Fi(x)| nowhere vanishes. f Hence the matrix Fi(x) has the desired 
properties. 

Let us now return to the proof of the theorem. Suppose that the matrices 
Qi> • • • , Qm have been so selected that the matrices Si, • ■ • , Sm-i determined 
by this choice have their elements in absolute value less than assigned positive 
quantities ai, • • • , am-i respectively within given circles of radii less than 
1 2^2 1 , • • • , \xm\ . We shall now show that Qm+i can then be so selected that 
the corresponding matrix Sm has its elements less in absolute value than a 
prescribed positive quantity am within any given circle of radius less than 

I "^m+l I • 

*Acta Mathematica, vol. 6 (1884), p. 229. 
t Sehlesinger, loe. cit., p. 23. 
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To this end we observe that the matrix 

Qm{x)VZllix) 

will not be singular except at x = Xm+i ■ Thus, by the lemma above, a 
matrix Qm+i{x) without finite singularities may be chosen so that the dif- 
ference matrix 

Qm+i ( X ) - Qmix) VzU ( X ) 

has elements of as small absolute value as desired for | x | ^ dm < \ Xm+i \ ■ 
Hence the same is true of 

Sm{x) = Qm+l{x)Vm+l{x)QZ^ (X) - I , 

as we see by multiplying the preceding expression on the right by Fm+i(x)Q,l' (.t). 

Thus we have a method of assigning Qm+i ( x ) as soon as Qmix) has been 
chosen, in such a way that the elements of Sm are less than a positive constant 
ttm f or I X I < (^m < I Xm+i | , whcrc ai + a2 + • • • is a convergent series. More- 
over, since lim Xm = » we may take lim dm = '^ also. 

But the elements of the infinite sequence of matrix products 

/4+i(x)Qr'(x) = [I + Smix)][I + Sm-iix)] ■■■ [Z + Si(x)] 

(.ni = 1,2, ■■■) 

will certainly converge absolutely and uniformly to the elements of a limit 
matrix if the elements of Smix) are less in absolute value than the mth term 
of a convergent series of positive terms.* Furthermore the determinant of 
this limit matrix is the limit of the product of the determinants of the factors 
and the sum of the logarithms of these factors will form a convergent series. 
Hence we infer that lim Rmix) is a matrix of entire functions E{x) of the 
form (5) and of determinant not zero save at Xi , .T2 , • • • at which points E (x) 
is equivalent to Vi{x) , F2 ( x ) Fi ( x ) , • • • respectively, so that E (x) has 
the desired properties. 

The most general matrix E{x) of entire functions with its singularities at 
Xi , X2 , • • • and equivalent to Fi ( x ) , F2 ( x ) Fi ( x ) , • • • at these points respec- 
tively has the form Q(x) E{x) tvhere Q (x) is a matrix of entire functions with 
no finite singular point. This fact is immediately evident. 

In the case re = 1 the matrix Qmix) , as determined by the method of the 
lemma, is merely an exponential with an exponent which is polynomial in x . 
Thus we have a thoroughgoing generalization of the factorization theorem in 
the case n = 1 , but one which appears much less explicit. Furthermore, the 
determination of each convergence factor Qm+i ( x ) is made to depend on the 
previous factor Qm{x) . 

* The proof for ra > 1 is not essentially more difficult than in the familiar ease n = 1. 
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4. Factorization of matrices of entire functions. Second method 

A more explicit form of factorization can be given by the aid of an extended 
elementary matrix 



(7) 



U(x) = 



where 

(8) 







1 •• 


• 


■«■(!)■' 




• 


1 •• 




I^'(^.)i 




■ 


•• 


. gP.(mo) fi _ ^\ . 


• 



Hm 



Prix) = X +-^ + "3 ■■■ +7' 



This matrix differs from the unit matrix only in some Ith column and reduces 
essentially to an elementary matrix for r = (compare with (1)). 

It is apparent that for a given value of r the lemma of § 1 may be extended 
to this type of matrix; i. e., a matrix F{x) with singularity of analytic type 
at X = a;o of multiplicity k may be written as Q{x) U (x) where Q has its 
elements analytic at Xo , and has a singularity of analytic type of order k — 1 
at cc = Xo if A; > 1 , or is not singular at .Tq if ^* = 1 . 

In fact, if we multiply F (x) on the right by a matrix 



Dix) 



1 • 


• <k) ■ 


•■ 


1 ■ 


■ <0 ■ 


• 


■ 


1 


• 







'"(^y 



which takes the same form as Z) of § 1 for x = .To and in which the Ith column 
contains the d's, we may reduce the elements of the Zth column of F (x) to 
zero at X = Xq (for some I). Consequently if we write 
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F{x)D(x) = Q{x) 



■ 1 


•■• 





1 ••• 





••• 


. 


••• 








pXa:/xo) 



\ xo) 












1 



the elements of Q(x) have the desired properties, since Q{x) is derived from 
F {x)D (x) merely by dividing the Zth column by 



„p,.(xlzot 



(-s> 



Multiplying through the above equation by the inverse of D (x) we have 
F (x) expressed as the product of Qix) and an extended elementary matrix 
U (x) like (7) with Ci = — di, Ci = — d^, • • • , c„ = — dn- 

It is to be especially observed that the choice of the quantities di and ac- 
cordingly d is independent of the choice of r which does not enter explicitly 
into D (x) for x = xq. 

We are now prepared to demonstrate immediately the following second 
factorization theorem for matrices of entire functions: 

Theorem III. Any matrix of entire functions E{x) with its singularities at 
the distinct points Xi , x^, • • • may be uritten in the form 

Eo{x) lim Um{x) Um-l (x) ••■ Ui{x) 

m=oo 

where E(,{x) is a matrix of entire functions of determinant nowhere zero and where 
Ui{x) , Uiix) , ■ ■ • are extended elementary matrices of the form (7) such that 
the finite singtdar point of U\{x) , Uiix) , ■ ■ ■ , Uk{x) lies at Xi , that of 
Uk+i (x) , • • • Ui{x) lies at x^ , and so on. 

Proof. The possibility of selecting extended elementary matrices Ui{x) , 
• ■ ■ , Uk{x) with singular point at .ti , such that the product Uk{x) ■ ■ • Ui{x) 
is equivalent to E{x) at x = xi , is immediately obvious. In fact we maj' 
determine the product by successive application of the extended form of the 
lemma of § 1, so that it is equivalent to E{x) at Xi. Consider now the 
product 

E{x)[lh{x) ■■■ L^(.x)r, 

which is a matrix of entire functions with its finite singular points at .T3 , 
Xz, • • • . Treating this matrix at x^ as we did E{x) at xi , we can find a second 
product of the form JJ]t+i{x) •■ ■ Ui{x) to which the product matrix under 
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consideration is equivalent at Xi . Proceeding in this way we obtain a sequence 
of matrices JJiix) , Uiix) , • • • . 

We can, in addition, secure the uniform convergence of the elements of 
the product matrix 

Um{x) XJm-i {x) ••• Ui{x) , 

for lim m = CO , within any arbitrary circle | a; | = p in the complex plane. 

For we have noted that, at each step in the determination of the successive 
factors, the constants Ci, Ci, • • • , Cn are chosen independently of the integer 
r . Thus, at each step after these constants have been chosen, we can take r 
so large that the elements of the Zth column of that factor, 

differ from the corresponding elements of the set 

0, 0, •••, 1, •••0 

within any assigned circle |x| = p <|xi:| by an arbitrarily small amount. 
Hence the elements of the corresponding factor Um{x) can be made to differ 
from the elements of the unit matrix by arbitrarily small quantities throughout 
this circle. 

It is therefore apparent that the factors Ui{x) , U^ix) , • ■ • can be so taken 
that the limit of their product exists uniformly within every fixed circle with 
center at the origin in the complex plane. The elements of this limit matrix 
Ei{x) are of course entire functions, and it has the desired equivalence proper- 
ties at xi , X2 , • • • . The matrix Eo{x) defined by the equation 

E{x) = Eo{x)Ei{x) 

is a matrix not singular in the finite plane. 

5. Equivalence problem for matrices of entire functions 

The results of § 3 permit us to solve the following problem: 
Equivalence problem. Given polynomial matrices Vi{x) , F2 ( x ) , • • • 
each with its only finite singular point at the distinct points Xi , Xi, • • ■ re- 
spectively; to construct a matrix E{x) of entire functions equivalent to Fi(x), 
F2 ( X ) , • • • at Xi, Xi, ■ ■ • respectively but with no other finite singular points. 

In fact, let us use the bracket notation [T {x)]^ as before to denote some 
polynomial matrix equivalent to7'(x)atx = ^ but otherwise not singular in 
the finite plane. Let us write 
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V^{x) = [V,(x)U = Fi(x), 

v,ix) = iv,{x)Vr'{x)u, 



Consider now the matrix of entire functions 

E(x) = \iinQm{x)Vm{x) ■■• Vi{x), 

jn=oo 

where the matrices Qm{x) are chosen as in Theorem II so that the limit is 
approached uniformly. It is readily seen that the product Vm(x) ■ ■ ■ Vi{x) 
is equivalent to F,- ( a; ) at x = Xi f or i = 1 , 2 , • • • , m. This is obvious for 
i = \ . Furthermore at x = X2 we may write 

F2(x) = A{x)V,{x)Vt'{x) 

where A{x) is not singular at x = xi . It follows at once that F2 ( x ) Fi ( x ) 
and hence F„ ( x ) • • • Fi ( x ) is equivalent to F2 ( x ) at x = X2 . Continuing 
in this way we may establish the truth of the statement made for i = 1 , 2 , 
• • • , m. 

The definition of E{x) now shows that it has the stated properties. Thus 
the equivalence problem has been solved. 

The most general solution of the equivalence ■problem is Eo{x) E{x) where 
Eo ( x) is an arbitrary matrix of entire functions of determinant not zero in the 
finite plane. 

6. Matrices with a finite number of essential singularities 

Thus far we have restricted attention to matrices with singular points of 
analytic type. We propose now to consider the question of the factorization 
of matrices with isolated essential singularities. We will restrict attention at 
first to the case in which there are only a finite number of isolated essentially 
singular points. 

I shall assume the truth of the following lemma: 

Lemma. Let L{x) he a matrix of functions for which x = 'Xi is an isolated 
essentially singular point. Then it is possible to write L{x) in the product form 
A(x) E{x) (x''' bij ) where A{x) is a matrix of functions not singular atx = <x> , 
where E{x) is a matrix of entire functions not singular in the finite plane, and 
where ki, k^, ■ • • , kn are integers. 

For the proof of this lemma the reader is referred to the Mathematische 
Annalen, vol. 74 (1913), pp. 122-133. With its aid we may at once prove 
the theorem which follows: 

Theorem IV. Any matrix F (x) with a finite number of isolated singularities 

Trans. Am. Math. Soc. 26 



400 G. D. bibkhoff: [July 

(essential or non-essential) at the distinct points X\, x^, • • • , Xm, x^+i = <» , and 
no others, can he written as a product of factors E{x) Wm {x) Wm^i (x) • • • 
Wx {x) where these factors W are 

Wi{x) = Ei(~-\{{x-xiY?bii]. 

\ X Xi J 

Here E (y) , Em (y) , • • • , Ei(y) are matrices of entire functions not singular 
in the finite plane and UP {I = 1,2, ■ ■ • , m; j = 1,2, ■ ■ ■ , n) are integers. 
Proof. It is clear that F [(xi + 1/^)] has an isolated essential singularity 
at I = 00 and so, by the above lemma, may be written as 

^i(^)-Ei($)(#5,,), 

where -Ei ( |) is a matrix of entire functions which is not singular in the finite 
plane. Transferring back to the variable x , we obtain the decomposition 

F{x) =F,{x)EA-^\[{x-x,r''Thii], 

where F\{x) is not singular at a; = a;i . It follows that Fi{x) has singulari- 
ties only at a;2 , • • • , Xm, and Xm+i = » . 
By the same process we find 

Fi{x) = F,{x)eA—^\[{x - X,)-''? hi,], 

\X — X2 / 

where Fi{x) is not singular save at 0:3, • • ■ , Xm, and Xm+i = «= • 

Proceeding successively in this manner we obtain a set of similar equations 
by introducing matrices Fz, • ■ • , F^. The last equation is 

Fm-i (x) = F,„{x)E,n{ — — r ) [ (a; - x™)-4""6iy] , 

\ X Xtn J 

where Fm{x) is not singular save at infinity, and hence is a matrix E{x) of 
entire functions of the desired type. 

From this succession of equations we get the stated decomposition at once. 
It is not difficult to infer that the above decomposition is essentially unique. 

7. The corresponding equivalence problem 

We can now pass at once to the solution of the corresponding problem : 
Equivalence problem.* Given matrices Gi{x) , • • • ,Gm{x) each with one 
isolated singularity at the distinct points xi, ■ • • ,Xm respectively in the finite plane 
and with no other singular points; to construct a matrix F (x) equivalent to 

*Cf. Proceedings of the American Academy of Arts and Sci- 
ences, vol. 49 (1913), pp. 521-568; in particular pp. 540-541. 



1916] INFINITE PRODUCTS OF ANALYTIC MATRICES 401 

matrices Gi{x) , • • ■ , Om{x) at x\, • • • , Xm respectively, and otherwise not 
singular save for a possible pole at infinity. 

Proof. We follow the method used in §§ 2, 5 in solving the similar problem 
for matrices of polynomials and entire functions. 

Suppose that T (x) has an isolated singularity at x = f . Let [r(x)]j be 
used to designate the matrix 



^(^)[(x -«)-*' 6.7], 



to which, in accordance with the lemma of the preceding section, F (x) is 
equivalent at x = ^ . This matrix has no other singularity save for a possible 
pole at X = CO . Here E(y) is an entire function of y and not singular in the 
finite plane. 

If now we write 

Hi(x) = [Gi(x)U = Gi(x), 

H^ix) = [G2{x)Hr'{x)U, 



H^{x) = [G^{x)Ht'(x) ■■■H-l,{x)l^, 

we infer at once that the matrix F (x) defined by 

Fix) = Hm{x) ■■■ Hi{x) 

has the stated properties (see § 2). 

The most general solution of the equivalence problem is Q{x)F (x) , where 
Q{x) is any polynomial matrix of constant determinant. 

8. The matrix with general singularities 

The factorization theorems of §§ 1, 3, admit of a further extension as follows : 
Theorem V. Let Xi , X2 , • ■ • be an isolated set of points in the finite complex 
p ane with a set of limit points [c], one of which lies at infinity in the plane. 
Further, let Ei{y) , Ei^y) , ■ • • be arbitrary matrices of functions of the form 
E{y){y''' 5ij ) where E{y) is a matrix of entire functions not singular in the 
finite plane, and ki, • • • , kn are integers. Then there exists a matrix G (x) of 
analytic functions given by a formula 



G{x) =Yimq,n{x)Ej-^\ 

m=ao \ X Xm / 



where Qm{x) is not singular in the finite plane, such that G{x) is equivalent to 
Em[l/{x — Xi)] • • • Ei[l/{x — X2)] at Xmim = 1,2, ■ • • ) and is otherwise 
not singular save at the points of [c] . 
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Proof. We shall follow closely the method used in proving the analogous 
Theorem II (§3). 

For convergence of a formula for G{x) of the type indicated, it is sufficient 
to establish that Qm+i (x) can be so chosen that, in a region outside of an 
arbitrary small distance of Xm+i and within a circle of arbitrarily large radius 
with center at this point in the complex plane, the elements of 



\ X Xfn J 



(9) Q™+i ( a; ) i?™+i :^-_— Q,l' ( x ) 

can be made to differ in absolute value from the corresponding elements of the 
unit matrix by less than an arbitrarily assigned positive quantity. 
But this is readily done as follows. The matrix 



^-^^^^M^^.) 



is obviously not singular except at a;m+i and at infinity. Hence by the lemma 
of § 3 we can find a matrix Qm+iix) , not singular in the finite plane, whose 
elements differ from those of the above matrix by less than an arbitrarily 
assigned quantity for 

r ^\X - Xm+l\ = R, 

r, R being arbitrary positive quantities. 

Consequently we may make the elements of the matrix (9) differ from those 
of the unit matrix by less than an arbitrary positive quantity in this region 
for r as small and R as large as desired. 

It is therefore clear that Qi(x) , Q2{x) , • ■ • may be taken in succession so 
that the limit in the theorem exists uniformly so long as x lies within any circle 
of arbitrarily large radius with center at the origin and not within an arbitrarily 
small distance of any point of [c] . 

The matrix G{x) thus defined has clearly the desired property. On the 
basis of this theorem we are led at once (as in § 5) to a solution of the problem: 

Equivalence problem. Given an infinite number of points xi,X2, ••• in the 
finite plane with a set of limit points [c], one of which lies at infinity, and matrices 
Ei{y) , E'i{y) , • • • of the form E{y) {y''' dij ) where E{y) is a matrix of entire 
functions, not singular in the finite plane, and ki, • • • , kn are integers; to con- 
struct a matrix G(x) equivalent to Ei[l/{x — xi)], E2[l/{x — Xi)], ■ • • at 
Xi,Xi, • • ■ respectively, and otherwise not singular save at the points of [c]. 

The most general solution of this equivalence problem is M (x)G (x) where 
M (x) is an arbitrary matrix of analytic functions whose singularities all fall at 
the points of [c]. 
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9. An APPLICATION 

I shall give now a simple application of the preceding theorems. This 
application may be stated as follows: Let e{x) be an entire function and let 
eo{x) , ei{x) , 62 {x) , ez{x) , ■•• denote the functions x, e{x), e{e{x)\, 
e{e[e{x)]] , ••• respectively. Furthermore let F (x) denote an arbitrary 
matrix of single-valued analytic functions. A matrix G{x) of single-valued 
analytic functions can then be found whose singular points fall at those of F[eic{x)] 
(k = 0, 1,2, •••) and at the limit points [c] of this set, and such that G{x) 
is equivalent to G [e{x)] F (x) save at these limit points. 

Let us write down for reference the symbolic infinite product to the left 

■■■ F[e2{x)]F[er{x)]F[eo{x)]. 

The singular points of any factor F [e„(a;)] not belonging to [c] are isolated 
singular points. 

Consider now any such singular point x not a point of [ c ] . It is possible to 
take a positive integer k so large that all the factors before the kth in the 
above product are not singular at x nor in its neighborhood, since x is not a 
point of [c]. 

We shall demand of a matrix G (x) that it be equivalent to the last k factors 
of the product 

A [ek-i{x)] ■■■A [eo{x)] 
at every point x . 

It is clear that if k be increased G{x) will still have the same property since 
the factors introduced are not singular at x, and appear on the left of the 
previous factors. 

The set of such points x is clearly numerable. In fact any set of isolated 
points is numerable. 

According to § 8 a matrix G{x) may be found which satisfies these relations 
of equivalence at the singular points x oi F [em{x)] {m = 0,1, •••) and which 
has no further singular points save the points of [c] . 

This matrix G {x) will have the desired property that G{x) is equivalent 
to G {e{x)] F (x) save at the points of [c] . 

In proving this statement it is first to be observed that if x is a singular 
point oi F [em(x)], and if a value y exists such that x = e{y) , then ^ is a 
singular point of F [Cm+i (j/) ] • For we have identically 

F[e„,ix)] = F[em+i{y)], 

where x = e{y) . Hence, if the matrix on the right is not singular at y , the 
elements of the equal matrix on the left are analytic in a; at x , since they can 
not be multiple- valued functions of x. Moreover, since as x approaches x, 
the variable y approaches y, the determinant of the right-hand side will not 
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approach zero in this case. Thus F [em{x)] would not be singular at x , 
contrary to hypothesis. 

The converse statement is obviously true, i. e., if F [ em+i ( 2/ ) ] is singular 
at y and if x = e{y) , then F [em{x)] is singular at x . 

Thus we derive the precise totality of singular points of successive factors 
of the symbolic infinite product from the singularities xo of the first factor 
F (x) by forming e_i ( xq ) , e_2 ( xo ) , • • ■ where e_i ( a; ) , e_2 ( a; ) , • • • denote 
the inverse functions to ei ( a; ) , 62 ( a; ) , ■ ■ • respectively, which are of course in 
general multiple valued. 

Consider now any point not belonging to [c]. 

If this point ^ is not a singular point of some matrix F [em{x)] then G {x) 
is certainly not singular at {, as has been stated. Moreover e({) is not a 
singular point for G (x) , inasmuch as in the contrary case some F [em{x)] 
is singular at e({) and thus the corresponding F [em+i{x)] is singular at f . 
Consequently G [e{x)]is not singular at ^ . It follows that G(x)is equivalent 
to G [e{x)] F (x) nt any such point { . 

It remains to consider a point x for which some F [em(.x)] is singular. At 
such a point G{x) has been chosen to be equivalent to a product 

F[e,{x)]F[e,-i{x)] ■■• F[eo(x)], 

where k has been taken so large as to include all the factors singular at x. 
The point e{x) will be a singular point for F [ e^-i ( a; ) ] unless m is zero. Put- 
ting aside this contingency for the moment, we see that G{x) is equivalent to 

F [e,_i{x)] F [e,_2{x)] ■ ■ ■ F leoix)] 

at e{x) . Consequently G [e(x)] is equivalent at x to 

F[e,{x)]F[ek-i{x)] ■■■ F[ei{x)]. 

Hence G (x) is equivalent to G [e{x)] F {x) at such a point x. 

If, on the other hand, m is zero then the integer k may be taken zero. In 
this case G (x) is equivalent to F (x) at x, while G{x) is not singular at 
e{x) . In this case also G{x) is equivalent to G[e{x)]F (x) at x. 

The italicized statement made is now proved. 

The usual existence theorems for linear difference and g-difference equations 
show that if e ( a; ) = a; + 1 or qx respectively, and if i^ ( a; ) is a polynomial 
matrix we may choose G{x) equal to G[e{x)]F {x) . 

Harvard University 
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